Superfluid to Bose-Glass Transition in a ID Weakly Interacting Bose Gas 
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We study the one-dimensional Bose gas in spatially correlated disorder at zero temperature, using 
an extended density-phase Bogoliubov method. We analyze in particular the decay of the one-body 
density matrix and the behaviour of the Bogoliubov excitations across the phase boundary. We 
observe that the transition to the Bose glass phase is marked by a power-law divergence of the 
density of states at low energy. A measure of the localization length displays a power-law energy 
dependence in both regions, with the exponent equal to —1 at the boundary. We draw the phase 
diagram of the superfluid-insulator transition in the limit of small interaction strength. 
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Interplay of disorder and many-body interactions is at 
the origin of a variety of interesting phenomena. As an 
example, a particle in a one-dimensional disorder poten- 
tial always displays localization [H . In a many-body sys- 
tem, on the other hand, derealization can arise as a con- 
sequence of interactions |2j. A related phenomenon in a 
repulsive Bose gas, is the quantum phase transition from 
a superfluid to an insulating Bose glass phase, as dis- 
order is increased @, [3, 0| ■ Trapped ultracold atomic 
gases are an ideal system for the investigation of the 
Bose-glass phase, thanks to the ability to tune the dis- 
order amplitude and the interaction strength [1, H, 0> H] • 
The achievement of the superfluid-Bose glass transition 
is however still under debate [tj [l(| • On the theoretical 
side, most efforts have been devoted to the strongly in- 
teracting regime of the phase diagram, through the study 



of the disordered Bose-Hubbard model [ll|, |12|, |13[ . Less 
attention has been paid to the weakly interacting dis- 
ordered Bose gas [3, EH EM E3, EH - a regime closer 
to that of recent experimental studies of Anderson lo- 
calization in an ultracold Bose gas 3, 2(j. The weakly 
interacting Bose gas in a continuous disorder potential is 
well described by a mean field approach, that provides 
closed- form expressions for the correlation functions. In 
low dimensional s yste ms, that are of particular experi- 
mental interest 0,13 IS] , this approach requires special 
care in the description of phase fluctuations, that play 
a dominant role by triggering the quantum phase tran- 
sition. Hence, a generalization of the theory to quasi- 
condensates (condensates with fluctuating phase) is nec- 
essary 
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In this Letter, we present a study of the ID disor- 
dered Bose gas within the extended number conserving 
Bogoliubov theory 2l| . The main physical quantities un- 
der investigation are the spatial coherence, the density of 
states (DOS) and the inverse participation number (IPN) 
of the elementary excitations. This latter is a good esti- 
mate of the spatial extent of the wave function, although 
it does not necessarily coincide with an exponential decay 
length of its tails We show that all three quantities 
allow to trace the mean field limit of the phase bound- 



ary between the insulating Bose glass and the superfluid 
quasi-condensed phases. We confirm the scaling of the 
IPN with energy, found by Gurarie et al. [IH and ex- 
tend it into the Bose glass phase. Surprisingly the gen- 
erally accepted fact 



131161 that the density of states 



is constant throughout the transition is contradicted by 
our numerical results: we find that the density of states 
of the Bogoliubov excitations diverges in the Bose glass 
phase. Furthermore we show that the loss of spatial co- 
herence in the Bose glass phase is dominated by weak 
links, across which the coherence drops sharply. Finally, 
we present numerical evidence that the critical disorder 
amplitude scales with the interaction energy as a power 
law. 

The N-body Hamiltonian describing the Bose system 
has the form 



H = dr 



&(r)H $(r) + |* f (r-)* T (r)*(r)*(r) , (1) 



where Ho = —h 2 df/(2m) + V(r) is the noninteract- 
ing Hamiltonian, ^(r) is the field operator, g is the 
coupling constant and V(r) is the disorder potential. 
Here we consider the case of a Gauss-distributed and 
Gauss-correlated random potential V(r), i.e. such that 



(V(r)V(r')) 



where A is the Gauss ampli- 



tude and i] is the spatial correlation length. 

In the mean field limit, it is useful to rewrite the annihi- 
lation operator in terms of a c-number quasi-condensate 
density po, and operators for the density and phase fluc- 
tuations, bp and 9 respectively: ^ = e %6 \J po + Sp. A 
correct definition of the phase operator, in the extended 
Bogoliubov method j2jj , requires the definition of a grid 
with step size I that fulfills the requirement on the to- 
tal density pi = (po + (Sp))l > 1. It can be shown that 
the quasi-condensate density obeys the Gross-Pitaevskii 
equation [24 1 



H + 9Po(r) V Po{r) = aVpoM, 



(2) 



where p, is the chemical potential. The density and phase 
fluctuations can be expressed in terms of the usual u and 
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v's of the standard Bogoliubov theory 25j, that obey the 
Bogoliubov-de Gennes equations 



Hq + 2gp {r) - M j Uj(r) + gp (r)vj(r) = E jUj (r), 

-9Po(r)uj(r) - (h + 2gp (r) - fi) vj(r) = E 3 v 3 (r). 

(3) 

The number-conserving formalism requires an orthog- 
onalization of the Bogoliubov modes Uj(r) and Vj(r), 
with respect to the quasi-condensate density po{r), 
thus obtaining the modes u±j(r) and v±j(r). Using 
Wick's theorem, the one-body density matrix G(r,r') — 
(&(r)4>(r')) takes the form 



G(r,r') = y/p(r)p(r')( 
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(r') 



(4) 

where at T = only the contribution from the quantum 
fluctuations appears. The approximations involved in de- 
riving Eq. Q require small density Sp/po <C 1 and phase 
fluctuations 59 /po <C 1. This mean field description 
holds in the high-density limit while in the low-density 
limit the regime of impenetrable bosons is reached and 
this approach is no more reliable. Throughout this work, 
we assume a spatially averaged quasi-condensate density 
par/ = 8 and, for the numerical calculations, we con- 
sider a finite system of length L and step size I. The 
simulations that follow are performed for L = 4096?7. 
We define the interaction energy U, through the rela- 
tion i7L|</>o(r)| 2 = gpo(r), with the normalized quasi- 
condensed wavefunction 4>q. In order for this discretized 
model to describe the continuous case, we always fulfill 
the basic requirement that the kinetic hopping energy 
t = be much larger than any other characteristic 

energy of the system. This must hold in particular for 
the energies U, A, p and E c = ■ This latter in turn 
implies that rj 3> I. In our simulations t/E c — 16, namely 
T] = 41. We solve the GPE using the Crank-Nicholson 
algorithm for the imaginary-time evolution of po- We 
then solve the linear Bogoliubov-de- Gennes problem by 
numerical diagonalization. Periodic boundary conditions 
are assumed both for the equations and for the randomly 
generated disorder potential V{r). 

We first study the long-range behaviour of the one- 
body density matrix. In the quasi-condensed phase this 
quantity is expected to have a power-law decay, as in the 
spatially uniform gas 21,23,[2y|. In presence of disorder, 
a transition to an exponential decay of the spatial corre- 
lation - characterizing the Bose glass phase - is expected 
below a critical value of U [3|. The quantity G(r, r') is 
affected by the specific shape of the disorder realization. 
Spatial average leads to the degree of coherence 

G(r,r') 
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FIG. 1: gi(r) in double logarithmic scale, computed for a sin- 
gle realization of the disorder at fixed A = 0.8E C , for different 
U. Inset: G(r,ro) in the Bose glass case (U = 0A8E C ). 



In our simulations, this quantity reveals to be self- 
averaging and reproduces directly the decay of the 
realization-averaged one-body density matrix. Our nu- 
merical analysis shows that the | ^ | 2 diverge as 1/Ej for 
Ej — > (in what follows, the zero of the energy scale is 
taken at the chemical potential p) . This behaviour seems 
to be generic, because it is also found in two limiting cases 
that allow for a simple analytical solution: the homoge- 
neous gas [24j and a Josephson junction in the limit of 
small tunneling (see below). Inspection of Eq. (j4]) 
then shows that the decay of long range correlations is 
driven by low energy excitations. To speed up calcula- 
tions, we therefore computed only the first N max = 2048 
eigenstates, from which we extracted the one-body den- 
sity matrix. We checked the convergence of the long 
range spatial coherence as a function of the cutoff N max 
and system size L. 

Fig. [T]shows the simulated gi(r) at varying U for fixed 
A = 0.8E C . The inset shows the quantity G(ro,r) com- 
puted for U — QA8E C , deep in the Bose-glass phase. We 
point out that our analysis of the long-range decay of 
gi(r) is carried out within the range [0, L/4], in order 
not to be affected by the periodic boundary conditions. 
The quantity gi (r) shows an exponential decay for small 
U, characterizing the Bose-glass insulating phase. By in- 
creasing U, it then increases at fixed r, up to a point 
where the interaction drives the system into a super- 
fluid phase, marked by a power-law decay. By further 
increasing U , the correlation at fixed r decreases again, 
in analogy with the decrease of coherence for increas- 
ing interactions in the spatially homogeneous case [24j ]. 
When increasing A at fixed U instead (not shown), the 
coherence always shows a monotonic decrease. 

The decay of correlations in the Bose glass phase is 
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FIG. 2: Averaged D(E) at fixed A = 0.8 E c for various inter- 
action energies. 



dominated by low-energy excitations having a phase flip 
character. Where a low energy excitation has a node, 
the coherence drops suddenly with a step that is related 
to the amplitude of this excitation. This behaviour of 
G(ro,r) is illustrated in the inset of Fig. [1] Physically, 
the two parts of the Bose gas on each side of the step can 
be seen as a weakly coupled Josephson junction. In the 
superfluid phase, no such abrupt drops in the coherence 
are found, but rather a smooth behaviour akin to the 
uniform system. 

It is interesting to note that the functional behaviour 
of the long range spatial coherence is determined by the 
interaction energy U and not by the density alone. In- 
deed, it is determined by the v's and the Bogoliubov-de 
Gennes equations (J3j) only depend on the product gpo. 
Eq. (fj| shows that the spatial coherence is reduced for 
decreasing density when keeping the interaction energy 
constant. 

We now turn to study the DOS of the Bogoliubov 
excitations, defined as D{E) = ^2jS(E — Ej). In the 
quasi-condensate phase, we expect D(E) to approach 
a constant for E — > 0, similarly to phonons in ran- 
dom elastic chains 12811 . For the Bose-glass phase, it 
has been argued 0, [l5| that the low-energy limit of the 
DOS should remain constant. Our results, however, do 
not support this behaviour in the mean field limit. We 
show in Fig. O the quantity D{E) plotted for 4 val- 
ues of U. For the largest U it clearly displays a con- 
stant limiting value for E — » 0, while D(E) develops a 
power-law divergence for the smallest U . The numerical 
accuracy of the data makes it difficult to extract the lim- 
iting behaviour of D(E) for intermediate values of U, al- 
though the power-law divergence in the Bose-glass phase 
is clearly assessed. We propose the following handwaving 
argument in support of the power-law divergence of the 



FIG. 3: Averaged L(E) for varying interaction at fixed A = 
0.8 E c - The horizontal line indicates a limit where the finite 
size effects start to play a role leading to saturation of L(E). 
Inset: Exponent of the power-law divergence of L(E) as a 
function of U/E c : the shaded zone marks the phase transition. 



DOS in the glass phase. By inspecting Eq.(TJJ for fixed 
r' = ro, we infer that the asymptotic behaviour at large 
r — ro is mainly determined by the term ■ |Ujx(?*)| 
in the exponent. Thus, an estimate of Eq.(f4|) in that 
limit is G(r, r ) ~ exp[/ \ve±(t)\ 2 D(E)dE] (see also 
Ref. (HI). Since ^.e-lO")! 2 always shows a l/E divergence 
for E — > 0, then the change of G(r, ro) from power-law to 
exponential must be determined by D(E) passing from 
constant to a power-law divergence. 

Another quantity of interest in connection with recent 
III 2(j[ is the localization length of the Bo- 
For a single particle, all 




experiments 
goliubov modes 
wave functions are expected to be exponentially local- 
ized in ID [29(, with the IPN monotonically increasing 
from a finite value at E — > — oo. In presence of inter- 
actions, as already pointed out, delocalized low-energy 
p/iase-excitations represent the major mechanism of long- 
range decoherence. The behaviour at E — > in the quasi- 
condensate phase should diverge as a power law E~ a but 

1. To 
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the value of a is still object of controversy 
study the localization of the excitations, we have com- 
puted their inverse participation number 



1 _ Jdr\v ±j (r)\ 4 
I, (fdr\v ±j (r)\* 



(6) 



and the corresponding realization-averaged quantity 

L ( E ) = Y.jWiE - E i)/ D { E )- We focus here 011 the 
Bogoliubov v± modes that are the only contribution to 

decoherence at T = 0. Our results for L(E) are dis- 
played in Fig. [3] Remarkably, we obtain a power-law 
divergence E~ a for E — ► 0, independently of the gas 
phase. The exponent a (plotted in the Inset of Fig. [3]) 
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varies continuously from a <C 1 for the lowest value of 
U, to a > 1 for the largest U considered. The finite size 
of the simulations limits the analysis for large U, deep 
in the quasi-condensate phase. This supports the sce- 
nario recently proposed by Gurarie et al. where a 
equals unity at the phase boundary and increases to 2 
for increasing interaction strength. We thus tentatively 
identify the a = 1 case as the phase boundary. More- 
over we find that a continuously decreases when going 
deeper in the insulator phase, apparently linearly van- 
ishing for U — > 0. This picture is consistent with the 
constant limiting value of L(E) at low energy, expected 
for the non-interacting case. 

The three quantities studied above, display a phase 
boundary for the same values of A and U, within the 
numerical accuracy. This allows to draw, in Fig. [4J 
a phase diagram for the quasi-condensate to Bose-glass 
transition at zero temperature close to the origin of the 
U — A plane. Symbols denote simulated systems in both 
phases. The shaded region contains the phase boundary 
and its width denotes the uncertainty in extracting an 
asymptotic behaviour from the finite system size. This 
uncertainty was reduced by the joined analysis of gi(r), 
the DOS and the IPN. In particular, it is remarkable 
how the boundary coincides with the a — 1 exponent 
for the IPN. The boundary appears to obey a power law 
A/E c = C(U/E c )i, with an exponent 7 = 0.75 ± 0.03. 
Finally, we remind that, at densities much lower than 
here considered, the mean field approximation becomes 
invalid. In the Bose glass phase, this occurs when the 
coherence is limited to a single maximum of the quasi- 
condensate density, the so-called Lifshitz glass phase 1J] . 
In the superfluid phase, the Bogoliubov approximation 
breaks down when the coherence is strongly reduced on 
the scale of the interparticle distance and an interaction 
dominated Bose glass phase is reached. The resulting 
reentrant behaviour of the glass phase as a function of 
interaction strength [i, 12 1 can consequently not be de- 
scribed within Bogoliubov theory. 

In conclusion, we have derived the mean field limit 
of the phase diagram for a ID Bose gas in presence of 
continuous, spatially correlated disorder. While analo- 
gous studies exist on the Bose-Hubbard lattice model, not 
much attention had been devoted to this limiting case. 
The first experimental works on the dilute Bose gas in a 



disorder potential ljj 20j essentially considered the case 
of vanishing interactions. In relation to the phase dia- 
gram, we have estimated for both these experiments the 
value U/E c ~ 1. In Ref. A/E c ~ 0.05. In Ref. 
the disorder parameter A/ J, corresponding to our A/f, 
is varied over a broad range. The continuous limit here 
considered is only appropriate to the situation A/t <C 1. 
The results here obtained provide useful indications for 
future experiments aimed at the characterization of the 
Bose glass phase boundary. 

We are grateful to I. Carusotto, D. Sarchi and P. Lugan 
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FIG. 4: Sketch of the phase diagram of the ID Bose gas as 
a function of interaction and disorder. (A): Bose glass; (□): 
quasi-condensate. The points denoted by a circle are those 
for which an asymptotic behaviour cannot be extracted from 
our simulations, implying proximity to the phase boundary. 
These points are therefore included in the shaded region that 
marks the transition. 



for enlightening discussions. This work was supported by 
the Swiss National Science Foundation through project 
No. 200021-117919. 



Electronic address: luca.fontanesi@epfl.ch 
1] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 
(1993). 

2] T. Giamarchi and H. J. Schulz, Phys. Rev. B 37, 325 
(1988). 

3] M. P. A. Fisher et. al, Phys. Rev. B 40, 546 (1989). 
4] R. T. Scalettar, G. G. Batrouni, and G. T. Zimanyi, 

Phys. Rev. Lett. 66, 3144 (1991). 
5] H. Moritz et. al, Phys. Rev. Lett. 91, 250402 (2003). 
6] L. Fallani et. al, Phys. Rev. Lett. 98, 130404 (2007). 
7] Y. P. Chen et. al, Phys. Rev. A 77, 033632 (2008). 
3] M. White et. al, Phys. Rev. Lett. 102, 055301 (2009). 
9] B. Damski et. al, Phys. Rev. Lett. 91, 080403 (2003). 
3] D. Delande and J. Zakrzewski, Phys. Rev. Lett. 102, 

085301 (2009). 
1] I. Bloch et. al, Rev. Mod. Phys. 80, 885 (2008); 

V. I. Yukalov, Laser Physics 19, 1 (2009). 
2] N. V. Prokof'ev and B. V. Svistunov, Phys. Rev. Lett. 

80, 4355 (1998). 
3] S. Rapsch, U. Schollwock, and W. Zwerger, Europhys. 

Lett. 46, 559 (1999). 
4] P. Lugan et. al, Phys. Rev. Lett. 98, 170403 (2007). 
5] P. Nisamaneephong, L. Zhang, and M. Ma, Phys. Rev. 
Lett. 71, 3830 (1993). 
[16] N. Bilas and N. Pavloff, Eur. Phys. J. D 40, 387 (2006). 



[1] 

[2] 

[3] 
[4] 

[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 

[12] 

[13] 

[14] 
[15] 



■5 



[17] P. Lugan et. al, Phys. Rev. Lett. 99, 180402 (2007). 
[18] V. Gurarie, G. Refael, and J. T. Chalker, Phys. Rev. 

Lett. 101, 170407 (2008). 
[19] J. Billy et. al, Nature 453, 891 (2008). 
[20] G. Roati et. al, Nature 453, 895 (2008). 
[21] C. Mora and Y. Castin, Phys. Rev. A 67, 053615 (2003). 
[22] J. O. Andersen, U. Al Khawaja, and H. T. C. Stoof, Phys. 

Rev. Lett. 88, 070407 (2002). 
[23] D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, 



Phys. Rev. Lett. 85, 3745 (2000). 

[24] L. Pitaevskii and S. Stringari, Bose- Einstein Condensa- 
tion (Clarendon Press, Oxford, 2003). 

[25] S. Shevchenko, Sov. J. Low Temp. Phys. 18, 223 (1992). 

[26] M. Schwartz, Phys. Rev. B 15, 1399 (1977). 

[27] G. S. Paraoanu et. al, J. Phys. B 34, 4689 (2001). 

[28] T. A. L. Ziman, Phys. Rev. Lett. 49, 337 (1982). 

[29] D. J. Thouless, Phys. Rep. 13, 93 (1974). 



